Accretion of nonminimally coupled scalar fields into black holes 
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By using a quasi-stationary approach, we consider the mass evolution of Schwarzschild black 
holes in the presence of a nonminimally coupled cosmological scalar field. The mass evolution 
equation is analytically solved for generic coupling, revealing a qualitatively distinct behavior from 
the minimal coupling case. In particular, for black hole masses smaller than a certain critical value, 
the accretion of the scalar field can lead to mass decreasing even if no phantom energy is involved. 
The physical validity of the adopted quasi-stationary approach and some implications of our result 
for the evolution of primordial and astrophysical black holes are discussed. More precisely, we argue 
that black hole observational data could be used to place constraints on the nonminimally coupled 
energy content of the universe. 
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I. INTRODUCTION 

The accretion of matter is one of the most studied 
physical process involving black holes. Assuming the va- 
lidity of certain energy conditions for the accreting mat- 
ter, the black hole mass will never decrease. In fact, if 
the null energy condition holds, no classical process can 
lead to mass decreasing for black holes The situation 
changes completely if quantum processes are allowed: a 
black hole can, in fact, shrink due to the emission of 
Hawking radiation Q . Such processes are particularly rel- 
evant, for instance, to Primordial black holes (PBH)[3]. 
One of the most striking features of PBH is that they 
could indeed evaporate completely due to the emission 
of Hawking radiation. It is known, in particular, that a 
PBH with mass smaller than the so called Hawking mass 
Mh = 10 15 g should have already evaporated by now. 
PBH with masses close to that limit are specially rele- 
vant because their emitted Hawking radiation might, in 
principle, produce observable effects in the present day 
universe [J]. 

The interest in these problems has increased con- 
siderably in the last years due to the many dark en- 
ergy phenomenological models that have been proposed 
to described the recent accelerated expansion of the 
universe Such models @ typically involve a scalar field 
pervading all the universe that could, in principle, be ab- 
sorbed by any black hole, implying consequently in new 
channels for black hole mass accretion Q- It is interest- 
ing to notice that the study of black holes growth in the 
presence of scalar fields has been initiated before [1] the 
discovery of the recent acceleration of the universe and, 
thus, before the proposal of any dark energy model. 

The mass evolution of any black hole is governed by 
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two competing processes. The first one is Hawking ra- 
diation, which decreases the black hole mass due to the 
emission of a thermal radiation. The other one, which 
tends to increase the black hole mass, is the accretion of 
the surrounding available matter and energy. The sur- 
vival or not of a PBH until nowadays, for instance, will 
depend on the detailed balance of these processes. The 
unexpected possibility that black hole masses could effec- 
tively decrease due to the accretion of exotic (phantom) 
dark energy@ was received with great interest because, 
mainly, it could alter qualitatively the evolution of any 
black hole, implying, occasionally, in observational con- 
sequences for both astrophysical and primordial black 
holes. Since phantom dark energy violates the usual en- 
ergy conditions, there is no contradiction between these 
results and the classical theory of black holes. Neverthe- 
less, one should keep in mind that the physical viability 
of models involving phantom energy has been constantly 
challenged by their severe inherent classical and quantum 
instabilities [lOj. 

In this paper, we study the mass evolution of 
Schwarzschild black holes in the presence of a nonmini- 
mally coupled scalar field. A quasi-stationary approach 
is introduced and the mass evolution equation is analyt- 
ically solved for generic coupling. Our main conclusion 
is that, for black hole initial masses smaller than a cer- 
tain critical value, one could indeed have mass decreasing 
even in the absence of the Hawking evaporation mech- 
anism and without any component of phantom energy 
in the model. This is a more robust scenario for mass 
decreasing of black holes due to the accretion of exotic 
matter since it is not plagued by the phantom energy in- 
stabilities. Moreover, one could have, in principle, mass 
decreasing for considerably larger black holes than the 
minimally coupled case, with possible implications for 
primordial and astrophysical black holes, which could be 
explored in order to place observation constraints on the 
nonminimally coupled energy content of the universe. 
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II. NONMINIMALLY COUPLED SCALAR 
FIELDS AROUND BLACK HOLES 

We are concerned here with a scalar field <j> governed 
by the action 



S 



d 4 x^[F( ( f>)R-d a ^d a ( f>-2V^)}, (1) 



surrounding a Schwarzschild black hole. Nonminimally 
coupled cosmological models of the ty pe (HI) have been 
intensively used in modern cosmology [llj. Models for 
which it is indeed possible to reach F(<p) = are known 
to be plagued with singularities The hypersurface 
F{4>) = marks, in a sense, the boundary between stan- 
dard (F(4>) > 0) and phantom-like (F(<fi) < 0) behavior 
for the scalar field c£[13j]. We are mainly interested here 
in models such that F(cf>) > everywhere since, in such 
cases, phantom-like behavior is excluded by construction. 

Since Schwarzschild spacetime is Ricci-flat, the equa- 
tion of motion for <f> obtained from ([1]) reads simply 



n<j> = v'(cf>), 



(2) 



and the associated energy momentum tensor is given by 



9ab 



(d c 0O c + 2V) + V a V fc F - g ab DF. 

(3) 

Note that, due to Ricci-flatness of Schwarzschild space- 
time, we have V;,T a b — 0. By adopting the usual 
Schwarzschild coordinates, the spherically symmetrical 
version of Eq. ([2]) will be given by 



d 2 4> 
dt 2 
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2M 



d_ 
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r J dr 
2M\ 



r J 



V'tt). (4) 



The standard formulation of the stationary Bondi accre- 
tion process [13] for this problem consists in considering 
solutions of (U) with the following boundary condition 



lim (j>{t,r) = (j> c {t), 



(5) 



where <f> c {t) corresponds to the cosmological homoge- 
neous and isotropic solution of the model ((T|), with cos- 
mological and Schwarzschild time coordinates identified. 
Since no back reaction of the scalar field is taken into ac- 
count, our approach requires that the energy content of 
the scalar field must remain bounded and small around 
the black hole. Once we have a solution <j)(t, r) of (H]) with 
bounded energy and obeying the boundary condition (0), 
we assume that its energy flux on the black hole horizon 
is completely absorbed by the black hole, implying that 



dM 
~dt 



=2M 



r z T t r dfl. 



(6) 



This problem was solved, for F{tf>) — 1 and V(<f>) = 0, 
in [15| . In the Eddington-Finkelstein coordinates (v,r), 



with v = t + r + 2Mlog(r/2Af — 1) corresponding to 
incoming light geodesies, the pertinent solution corre- 
sponds to the stationary configuration 



v,r)=(3 + r y[v — r + 2M log 



(7) 



with (3 and 7 constant. We do not expect to have sta- 
tionary solutions like this for the generic model (fT]). In 
fact, stationary solutions are possible only for actions 
that are invariant under shifts <j> — > (f> + A, see [l6|. We 
can, however, adopt a quasi-stationary approach based 
on the observation [l?) that, for slowly varying cosmo- 
logical solutions <frc(t), the "delayed" field configuration 
given by 



(v, r) = 4> c 



2M log ■ 



2M 



(8) 



is an approximated solution of Q for certain potentials 
V((j)). The validity of this approximation will assure, of 
course, the validity of our quasi-stationary approach. By 
substituting §8§ in ^ one gets 



2M /2M\ 2 /2M x3 ' 



V r 



V r 



4> C +V'(<f> c ) = 0, (9) 



with the dot standing for the derivative with respect 
to t. Hence, our approximation is valid if <p c w and 
V'(<p c ) w 0. Due to the typical cosmological time scales, 
the assumption of a quasi-stationary {<j) c 0) evolution 
around the black hole is not, in fact, too restrictive. The 
same is true for the assumption V'(4> c ) w 0, but the argu- 
ment is more involved. Assuming a small variation of </> c , 
the potential can be linearized as V(<p c ) = /i</> c , since the 
constant factor is irrelevant here. In this case, equation 
will be a linear equation, and it is possible to find a 
stationary solution obeying the Bondi boundary condi- 
tion (J51). The approximation will be valid provided <f> c is 
small and r is kept smaller than the cosmological horizon 
scale, see Q for the details. It is interesting to notice that 
the explicit examples of failure of the approximation © 
presented in [l7j corresponds clearly to situations where 
one cannot assure <j> c or V'(<fi c ) » 0. 
For the solution (15]), one has 



\ r J 



AM' 



(10) 



Also from ([8j), we see that, on the black hole horizon, 
the field cj> assumes the value of <p c , propagated along 
a incoming light geodesic, but arriving with a certain 
"delay" 17 . Our quasi-stationary analysis neglects also 
such delay and, hence, in the quasi-stationary approxi- 
mation 

0c(*) «<?W + 0oo(*-*o), (11) 
with and 0oo constants, we have 

M = 16ttM 2 (1 + F") 4^ - 4ttMF'0 oo . (12) 
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For the minimal coupling case, F(<fi) = 1 and (fT2|) re- 
duces to the usual scalar field accretion rate[15||- It is 
clear, however, that for the nonminimally coupled case 
one could have, in principle, M < even in the absence 
of phantom modes. The rate ([12]) corresponds only to the 
accretion of the scalar field. The complete mass evolution 
equation is obtained by adding to the right-handed side 
a term cx M~ 2 corresponding to the Hawking radiation. 
As we will see in the next section, the fact that the two 
accretion terms in (|12p have different signs and differ- 
ent powers of M will imply in the existence of a critical 
mass M cr delimiting the mass increasing and decreasing 
accretion regimes. 

We finish this section by noticing that the possibility 
of negative energy fluxes for nonminimally coupled scalar 
fields and their implications for mass decreasing process 
involving black holes has been already considered previ- 
ously in another context, namely in the investigation of 
the generalized second law of thermodynamics [18j. 



the scalar field. In fact, in this case, the denominator of 
(TTSJ) vanishes for t — t cr , with H(t cr ) = Mq , implying 
that the black hole grows up to infinite mass in a finite 
time. The larger is the black hole initial mass Mo, the 
shorter is t CI . In contrast to the < Mo < M cr case, such 
behavior for Mq > M cr is similar to that one observed for 
the minimally coupled case F = 1. The qualitative evo- 
lution for the case Mo = M cr will depend on the details 
of the functions f(t) and g(t). 

For situations with large M cr , the inclusion of Hawking 
radiation will alter qualitatively only the final instants of 
the mass decreasing process. In such a case, for Mo < 
M cr , the black hole also disappears, but now in a finite 
time, since Hawking radiation dominates the process for 
M(t) <C 1. In fact, for M > M cr , the Hawking radiation 
term can be neglected and the dynamics are essentially 
that one described by (|15|) . Let us now consider some 
explicit examples of the coupling function F(<f>) in order 
to elucidate these points. 



III. MASS EVOLUTION 



A. F(</>) = !+£0 



For a generic coupling function F(<j)), the complete 
mass evolution equation has the general form 



M = f(t)M 2 - g(t)M 



a 



(13) 



where f(t) and g(t) are smooth functions and a is a char- 
acteristic constant for Hawking radiation. Let us con- 
sider, initially, only the accretion process (a = 0). By 
introducing M(t) = G(t)P(t), with 



G(t)=e-& BWdB , 



(14) 



we obtain a separable equation for P(t), which can be 
easily solved leading to the following solution for l|13p 
with a = 



M(t) = 

where M (to) = Mq and 



M Q G(t) 
1 - M H(t) ' 



H(t)= / f(s)G(s)ds. 



(15) 



(16) 



Typically, if the denominator of (|15[) does not vanish, 
the mass M(t) decreases according to (Tl4l for positive 
g(t). Mass increasing solutions appear when the denom- 
inator vanishes. For positive and well behaved f(t) and 
g(t), the function H(t) will be monotonically increasing 
and bounded by = limt_,oo H(t), leading to a criti- 
cal mass M cr = H^, 1 . Any black hole with initial mass 
Mo such that < Mo < M cr , even in the absence of 
Hawking radiation, will disappear due to the accretion 
of the scalar field, but such process typically will take an 
infinite amount of time. On the other hand, those black 
holes with initial masses Mq > M cr will grow by accreting 



In this linear coupling case, equation (|13p is au- 
tonomous, with f(t) — Wircp^ and g(t) = 47r£</> 00 , and 
can be integrated by quadrature for any value of a. We 
do not need, however, the exact solution here. We assume 
£ and <j> to be both positive in order to avoid possible 
singularities [l3| and, without loss of generality, <o = 0. 
The functions Git) and H it) are in this case 



Git) 



and 



Hit) 



40 c 



For (paa positive, we have 



M, 



(1 - Git)) 



and 



cr — ^Too i 



1 M 
■log- 



4tt£0oo Mq - M c 



(17) 



(18) 



(19) 



(20) 



Notice that, for typical cosmological situations, (fioo is 
small, implying in large values of M cr for £ of the order 
of unity (in Planck units). In these cases, the Hawking 
radiation is important only in the final instants of the 
mass decreasing phase. 



B. 



F(4>) = 1 + ^ 2 



We assume £ > 0. We have /(<) = 16tt(1 + 2£)<^ c and 
g(t) = 87r£ (<poo<P<x> + 'fitotj in this case. The pertinent 



4 



functions are, for to = 0, 

G(t) = e^^C 2 ^ 00 ^ 00 *" 1 "^' 3 ) (21) 

and 

H(t) = 16 7 r(l + 2£)0 2 o / e^ 4Tr «( 2 ^^ s+ ^ s2 ) da. (22) 



The critical mass is given by M cr = H^ 1 , with 



= 4tt 



1 + 2$ 



1-o-erf 



(23) 

where <r = sgn^oo and erf(x) is the error function [l^]. 
For the typical cosmological situations we have that (poo 
is very small, leading to 



4tt(1 + 2C) 

Notice that, as in the previous case, M c[ oc 4>^ ■ 



(24) 



C. F{<t>)=e 



t<t> 



2g£(0oo+0oot) U2 



In this case, we have /(t) = 167r 6-+ £ 
and = 4 7 r£0 oo e^~ +< ^ t \ leading, for t = 0, to 

G{t) = exp (-47re^~ (e^°°* - l)) (25) 



and 



(t) = 16tt^ jf (l + fV w - + *-' ,) ) G( a )ds. (26) 



The critical mass is given by 



A-r 1 = 



+ exp (47re^°° ) T (0, 4ne^ 



(27) 

where T(z, x) is the incomplete Gamma function [191]. For 
<ioo small, we have 



a + 4£ 



2 ^oo ' 



(28) 



where a is a numerical constant of the order of unity, 
namely a = 167re 4l T(0, 47r) « 3.72. Again, we observe 



the same behavior M cr oc t^ 1 . 



D. Radiation era with F{cj>) = 1 + £0 

The previous examples involve only the nonminimally 
scalar field in the quasi-stationary approximation. This 
is not enough, for instance, to describe PBH, since they 
were created in the primordial universe and have existed 



for eras where dark energy was not the gravitationally 
dominant content of the universe. In the radiation domi- 
nated era, in particular, the universe was filled and dom- 
inated by ultra relativistic matter which energy density 
is described in Planck units by 



(29) 



327rf 2 ' 



Such an energy density has been also available to be ac- 
creted by the black hole and should be incorporated in 
our analysis. The case of linear coupling F(<p) = 1 + £0 
in the presence of radiation with energy density (f2"9")l cor- 
responds to the choices f(t) = Wncp^ + (3/2)i~ 2 and 
g(t) = 47r£0 oo . The G(t) and H(t) functions in this case 
are 



G(t) = e -4*«*«>(*-*o) 



(30) 



and 



H(t) = ^-(l-G(t)) + |jf 
leading to 



fl -a e -4 fl rW»(.-to) ds 



to 



with (3 = 47r£</> 0o t . Since 



lim xlY— 1, x) = 1, 
x— >o 



we have in the present case 



M cr = H-i » -L ( 1 ■ 3i 



40 o 



&4>octo 



(31) 



(32) 



(33) 



(34) 



if P is small. 



IV. DISCUSSION 

If we assume that (p^ is of the same order of the critical 
density of the universe today (po ~ 10 _29 g/cm 3 ), we 
have M cr » 10 56 g for coupling constants £ of the order of 
unity (in Planck units) in the three first cases considered 
in the last section, allowing all the black holes in the 
universe to be in the shrinking phase today. In fact, 
even the galactic supermassive black holes (SMBH) with 
M w 10 6 M Q m 10 39 g are far below such a limit. These 
black holes would be shrinking today according to (|15p . 
The exact characteristic decaying time will depend on the 
particular coupling function. For the case of the linear 
coupling, the characteristic time is, according to (| 1 T[) . 
10 17 s, similar to the universe age. Notice that all the 
other coupling functions considered in the last section 
lead, typically, to faster decreasing mass regimes. 
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The fact that there are likely many black hole around 
us might be used to constraint the nonminimally cou- 
pled energy content of the universe during the cosmo- 
logical history. Let us consider, for simplicity, the last 
example of the previous section: the linear coupling case 
during the radiation dominated era. Suppose that the 
dark energy content of the universe has changed lightly 
after, say, to — Is. In this case, 4>ooto ~ 10 -18 in Planck 
units, justifying to take (3 w in (|32|) and leading to 
M CI w 10 38 g for a coupling constant £ of the order of 
unity. Thus, only PBH with mass greater than 10 38 g 
would escape from the shrink phase. Notice that this 
mass is extremely large if compared with the usual Hawk- 
ing mass Mh = 10 15 g. Observational constraints on the 
PBH mass cutoff[4] could be used, in principle, to estab- 
lish constraints on the non-minimal coupling parameter 
£, although the details depends on the coupling func- 
tion F(tfi). If we take to = 10 11 s, corresponding to the 
radiation-matter equality era, we will have </>ooto * 10~ 7 , 
leading to M cr w 10 49 g. This is, again, a huge mass and 
implies that virtually all black holes present at the end 
of radiation era have existed during all the matter dom- 
inated era in a shrinking regime. They should have lost 
two thirds of their mass by now, suggesting that observa- 
tional data about SMBH could also be used to constraint 
the nonminimally coupled energy content of the universe. 

We finish by noticing two points. First, one knows that 
it is not expected, in general, to have constant values for 
4>oa and along the cosmological history. Equation 
(|13p accommodates also situations where ipoa and (j>oo 
are functions of t. However, we should keep in mind that 
our formalism is based on the assumption of a quasi- 
stationary evolution, requiring cj> c (t) ~ in order to work 
properly. One needs to take backreaction into account in 
order to treat non stationary situations, see, for instance, 



[201 ] for a recent discussion. 

The second point is related with the hypothesis that 
cj) is a field test around a Schwarzschild black-hole. This 
is a good approximation provided that the energy con- 
tent of the scalar field (dark energy) is negligible when 
compared with the black-hole Physics scale. For the 
much larger cosmological scale, on the other hand, the 
scalar field is indeed the dominant energy content, be- 
ing the sole responsible for the accelerated expansion of 
the universe, usually described by a quasi-de Sitter so- 
lution. In our universe, these two scales are very differ- 
ent. Since the dark energy content is so small, in order 
to probe the quasi-de Sitter properties of the spacetime 
one needs to consider length scales of the same order 
of the Hubble radius. It is perfectly possible, in par- 
ticular, to apply condition Q in a region far from de 
black-hole (large r), but still far from the cosmological 
horizon. Furthermore, provided that the effective cos- 
mological constant of the accelerated expansion is small, 
the dynamics near the black-hole horizon are essentially 
the same of the Schwarzschild case, implying that (TT!?|) is 
still valid. From a theoretical point of view, however, it is 
certainly interesting to consider the problem of accretion 
onto Schwarzschild-de Sitter black-holes as it is done, for 
instance, in [2l| for the case of perfect fluids and min- 
imally coupled fields. We already know, however, that 
our present analysis should arise naturally in the limit of 
small A. These points are now under investigation. 
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